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Abstract 

The critical endpoint (CEP) and the phase structure are studied in the Polyakov-loop extended 
Nambu-Jona-Lasinio model in which the scalar type eight-quark (cr^) interaction and the vector 
type four-quark interaction are newly added. The interaction largely shifts the CEP toward 
higher temperature and lower chemical potential, while the vector type interaction does oppositely. 
At zero chemical potential, the a'^ interaction moves the pseudo-critical temperature of the chiral 
phase transition to the vicinity of that of the deconfinement phase transition. 
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The position of the critical endpoint (CEP) in the phase diagram is one of the most 
interesting subjects in hot and dense Quantum Chromodynamics (QCD). With the aid of 
the progress in computer power, lattice QCD simulations have become feasible for thermal 
systems at zero chemical potential (/i) [l]. For finite chemical potential, however, lattice 
QCD has the well-known sign problem, so that only a few works were made to determine 
the position of the CEP j^, ^. 

As an approach complementary to first-principle lattice simulations, one can consider 
several effective models. One of them is the Nambu-Jona-Lasinio (NJL) model ^. In the 
original NJL model that includes scalar and pseudo-scalar type four-quark interactions, it 
was found that there exists a CEP in the phase diagram 5, ^. However, the CEP is located 
at a lower temperature (T) and a higher /i compared with the one predicted by a lattice QCD 



simulation [2I] and by the QCD-like theory 







Moreover, recent empirical analysis |9|] of 



T]/s, the ratio of shear viscosity to entropy density, suggest there is a CEP at Tg ~ 165 MeV 
and /ie ~ 50 — 60 MeV that is much higher T and lower /x than the NJL model predictions. 

Kashiwa et al. lOj showed that in the NJL model the scalar type eight-quark (cr^) inter- 
action newly added shifts the CEP toward higher T and lower /i. However, the location is 
still far from predictions of lattice QCD, the QCD-like theory and the empirical analyses. 

It is also known in the chiral hadron model 
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12l | that a CEP appears at T much 



higher than the prediction of the NJL model. This implies that the deconfinement phase 
transition plays an important role in determining the position of the CEP. Although the 
NJL model is a useful method for understanding the chiral symmetry breaking, this model 
does not possess a confinement mechanism. As a reliable model that can treat both the 
chiral and the deconfinement phase transitions, we can consider the Polyakov-loop extended 
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14, 



15 



16|, 



17 



y 



19|, 



21 



22]. In the PNJL model the 



NJL (PNJL) model 

confinement /deconfinement phase transition is described by the Polyakov loop. Effects of 
the Polyakov loop make the CEP move to higher T and lower /i than the NJL model 
predicts [3] • The position of the CEP is still far from the predictions of lattice QCD and 
the empirical analyses. 



Meanwhile, .t wa= recently reported that the vecto. type interaction Hfl ESl i. 

necessary to realize the heavy neutron star [26'|. This may indicate that the vector type 
interaction is necessary in the finite n region of the phase diagram. 



In this letter, we study effects of the and the vector type interactions on the position 
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of the CEP and the interplay between the chiral and the deconfinement phase transitions, 
by using the PNJL model with these interactions. 

The model we consider here is the following two-flavor PNJL model with the vector type 



interaction 
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24 



25l | and the cr^ interaction 
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281,1291,13 



£ = q{t^^D'^-mo)q + C,nt-KmA],^[A\,T), 



(1) 



where q denotes the two-flavor quark field, mo the current quark mass and D'^ = — iA^ 
the covariant derivative. The field A'^ is defined as A'^ = 5^ogA'^^ with the gauge field A^, 
the Gell-Mann matrix and the gauge coupling g. The interaction of the NJL sector, C^^t, 
is 



(2) 



where f stands for the isospin matrix, and Gs, and Gss denote coupling constants of 
the scalar type four-quark, the vector type four-quark and the o"^ interactions, respectively. 
The Polyakov potential U, defined in Eq. ( ITSl) . is a function of the Polyakov loop ^ and its 
conjugate 3, 



(TrL), <P 



(TrLt), 



(3) 



with 



L(x) = Pexp i / (iryl4(x, r) 



(4) 

where V is the path ordering and A4 = iAo- In the chiral limit (mo = 0), the Lagrangian 
density has exact SU{2)i^ x SU{2)^ x f/(l)v x SU{3)c symmetry. 

The temporal component of the gauge field is diagonal in the flavor space, because the 
color and the flavor spaces are completely separated out in the present case. In the Polyakov 



gauge, the Polyakov loop matrix L can be written in a diagonal form in the color space 



— (TiL) = — fTr e^^^'^^'^^'^'^^'^^h 



— (Tr diag(e*"'"-', 

i Vp 



(5) 
(6) 



where 0a = + ^s/v^, (ph = -h + ^s/v^ and <^c = -(</>a + </>b) = -2(f)8/V3. The 
Polyakov loop is an exact order parameter of the spontaneous 1j{Nc) symmetry breaking 
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in the pure gauge theory. Although 'L{Nc) is not an exact symmetry in the system with 
dynamical quarks, it still seems to be a good indicator of the deconfinement phase transition. 
Therefore, we use ^ to define the deconfinement phase transition. 

Under the mean field approximation (MFA), the Lagrangian density becomes 



(7) 



where 



(7 

u 



{qq), p^{T,fi,a,^,3) = (qjoq) , 



(8) 
(9) 
(10) 



In the 1/Nc expansion, the eight-quark interaction after the MFA, qFq^qF'q)^ is of order 
and accordingly the same order as the four-quark interaction after the MFA, qrq{qr'q) , 
where F and F' are vertex matrices. Therefore, we can not ignore higher multi-quark 
interactions in general. As a starting point, we take into account the scalar type eight- 



quark (cr"^) interaction that is known to affect the position of the CEP [l^. Using the usual 



techniques 3l|, |32| , one can obtain the thermodynamical potential 



-T\nZ 



-2N,V 



(27r)3 



3E(p) + -{Tr,ln(l + Le-'^^"(p)) 



+Tr,ln(l + L^e~'^^^(p))} + |f/(a,cj) + W(^, T)|y, 



(11) 



(12) 



where E{p) = a/p^ + M^, = E{p) ± p ,M = rna + and p = p + E^. After some 
algebra, the thermodynamical potential i7 becomes 17 1 

d^p 



n = -2NfV 



(2ny 



+i In [1 + 3{<P + <Pe-'3^'(P))e-^^"(P) + e'^'^^^^P)] 
P 

+i In [1 + 3i§ + <Pe-^Enp)y-PEnp) + e-3PE+ip)^ 
+ {U + U)V. 



(13) 



In the T = limit, the PNJL model is reduced to the NJL model, since the Polyakov 
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loop is included only in U that has no a dependence, as shown in 



-E-{jp))E-{y>) 



(27r)3 

[/(a, pv(r ^ 0, /i, a)) + W(T ^ 0, ^) 



(14) 



We use U of Ref . 
at finite T 



331, 



34: 



17| that is fitted to a lattice QCD simulation in the pure gauge theory 



U_ 

2^4 



2 6 ^ ^ 4 ^ ^ ' 



&2(T) = ao + ai(^^) +a2(^Y) +"3( 



(15) 
(16) 



where parameters are summarized in Table I. 

The Polyakov potential yields a deconfinement phase transition at T = Tq in the pure 
gauge theory. Hence, Tq is taken to be 270 MeV predicted by the pure gauge lattice QCD 
calculation. 



ao 


0.1 


a2 


^3 


^>3 


hi 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



TABLE I: Summary of the parameter set in the Polyakov sector used in Ref. 
are dimensionless. 



1^. 



All parameters 



Since the NJL model is nonrenormalizable, it is needed to introduce a cutoff in the 
momentum integration. In this study, we use the three-dimensional momentum cutoff 

rf3p 1 rA 



p'dp. (17) 



(27r)3 27r2 

Hence, the present model has five parameters mo. A, Gg, Gss, and in the NJL sector. We 
simply assume mo = 5.5 MeV. In the case without the interaction, we use A = 0.6315 
GeV and Gg = 5.498 GeV^^, which reproduce the empirical values of the pion decay constant 
and the pion mass, /^r = 93.3 MeV and M^^ = 138 MeV. In the case with the o"^ interaction, 
we take A = 0.6315 GeV, = 5.002 (5.276) GeV^^ and G.gt^o = 0-2476 GeV^^ (0.1109 
GeV~^), which reproduce the pion decay constant = 93.3 MeV, the pion mass = 138 
MeV and the sigma meson mass M„ = 600 MeV (650 MeV) 
larger Ggg yields smaller M^-. 
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It should be noted that 



For the vector coupling constant Gv, we take three values, 0, 0.25Gs and O.SGg. The case 
Gv = O.SGs is obtained by one gluon exchange calculation in perturbative QCD jsst. The 
case Gv =0.25Gs is obtained by the instanton-anti-instanton molecule model 24 1. 

Table II summarizes the parameter sets we take. 



model 


Gs 




Original 


5.498 





Original + cj^ 


5.002 


0.2476 


Original + cj^ {M^ = 650 MeV) 


5.276 


0.1109 



TABLE II: Summary of the parameter sets of the NJL part. The coupling constants are shown in 
GeV~2. For ah cases, we take A = 0.6315 GeV and ao = a{T = 0, fj, = 0) = -0.03023 GeV^. 



Stationary conditions for cr, p^, <P and <P become 
1 df2 



with 



V da 



1 dQ 



Gla + 6iVfG: X 



f d^p f M r h(-\a,<P,§;T,i2) h(+\a,<P,§;T, p) ^ 
J (27r)3tE(p)l gi-)ia,<P,3;T,p) gMia,<P,§;T, p) ) ' ^ ^ 



2G,p, + 12NiG, X 



(27r)3 lgi~){a,<P,<P;T,p) gi+){a,<P,<P;T, p) 

1 df] 



0, (19) 



V d$ 2 



r d^p f e-^^^'(p) e-^^^(p) 
~^^f^y (2^\^{-)(a,<?,^;r,/i) + ^(+)(a,^,^;T,/i) l " °' 

1 dil rp4. _ 

— (-62 (r)<?- 63^2 ^64^2^) 



V d<P 2 



d^p f e-^^-(P) e-^^^^(p) 
/ (2^l^{-)(a,^,^;T,/x) + ^^^(a, ^, ^; T, p) I " °' ^^^^ 



(7(-)(cr, ^, ^; T, /i) = 1 + 3{<P + ^e-^^"(P))e-^^"(P) + e'^'^^'^P^ (22) 

g^+\a,$,§;T,p) = 1 + 3(^ + ^e-^^^(P))e-^^''(P) + e'^^^^^^P), (23) 

h^'\a,^,3;T,p) = <Pe-^E'i^) ^^^^-2pE-i^) ^^-WE-(v)^ (24) 

/i(+)(a, ^, ^; T, p) = c?e-'^^''(P) + 2^e-2^^^(P) + e'=^^^^(P), (25) 
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where the effective couphng G* is defined as G* = 2Gs + 12Gs80'^. 

Values of a are directly determined from minima of the real part of the thermodynamical 
potential, while $ and Pv are obtained by solving Eqs. (ITS]) - (PT]) . 



Following Refs. 



14 



2l|, we define the susceptibilities as 



X 



( 



1 



detC 



^-1 



(26) 
(27) 



where C is the matrix of dimensionless curvatures 

/ /3 d'^Q P d'^Q 



c 



\ 



_J3d'^Q 

13 a^si 



' 2Gs/l^ dcrd-P 

p 



p a'^n 

'2GlA^daa<P 

p a'^Q 
A3 a<pa<p 

p a'^n 

A3 g<p2 



\ 



(28) 



2Gs/i2 a'paa A3 a'Pd'p 

In this study, the susceptibility is used to determine a pseudo-critical temperature of 
crossover. 



(a) 

Chiral Condensate 



PNJL model 
PNJL + o^model 




lO 
O 

X 



X 
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(b) 


PNJL model 




PNJL+o* model ■ 











0.25 
T [GeV] 



0.3 



3.2 



0.25 
T [GeV] 



0.3 



Fig. 1: The left panel graphs the T dependence of the chiral condensate {qq) / {qq)T=iJ,=o and the 
Polyakov loop at fi = 0. The right panel graphs the chiral and Polyakov loop susceptibilities. 

First we discuss the T dependence of the chiral condensate, the Polyakov loop and their 
susceptibilities in the case of = 0. Figure 1 shows two results of the original PNJL and 
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the rNJI+cr"^ models, since the vector type interaction does not contribute to the phase 
transition when fi = 0. From the left panel, we can see that the interaction makes 
the chiral phase transition sharper and its pseudo-critical temperature lower in the PNJL 



model as well as in the NJL model 
model 



Similar effects are also seen in the three flavor NJL 



29 



30|. Meanwhile, the cr^ interaction affects the Polyakov loop little. As a result. 



the pseudo-critical temperature of the chiral phase transition goes down to the vicinity of 
that of the deconfinement transition. 

The pseudo-critical temperatures are calculated also in Refs. 



20(1 but not as that of Ref. 



17 



20| . Our regularization 



171], that is, in the present 



Doth the vacuum and T-dependent terms in the 



17( 1 the cutoff is made only for the vacuum term. 



scheme is the same as that of Ref. 
work the momentum cutoff is taken for 
square bracket of Eq. f[T^ . while in Ref. 
Consequently, our result is consistent with that of Ref. [20], but somewhat deviates from 
that of Ref. (jj. 

The pseudo-critical temperatures (Tc) can be clearly defined by the peak of the suscepti- 
bilities, as shown in the right panel of Fig.l. In the original PNJL model the pseudo-critical 
temperatures of the chiral and the deconfinement phase transitions are Tc = 273 MeV 
and 242 MeV, respectively. In the PNJL model, the corresponding values are = 247 
MeV and 236 MeV. The differences between the chiral and the deconfinement pseudo-critical 
temperatures are 31 MeV in the original PNJL model and 11 MeV in the PNJL+a^ model. 
Thus, the cx^ interaction makes the difference smaller. 

Second we discuss the behavior of the chiral condensate and the Polyakov loop near the 
CEP (Te,/ie)- Figure 2 shows the T dependence of the chiral condensate, the Polyakov loop 
and their susceptibilities near the CEP. As clearly seen in the right panel, both the chiral 
and Polyakov loop susceptibilities diverge at the CEP. This indicates that the two phase 
transitions are second order at the CEP. Further discussions on the CEP will be made in a 
forthcoming paper. 

Finally we show the phase diagram in Fig. 3. In the region of T < Te and /i > /le the 
chiral and deconfinement transitions are first order and occur at the same time. This can be 
understood by the generalized Clausius-Clapeyron relation for systems with multiple order 



parameters, which ensures that t 
the transitions are first order 



leir discontinuities appear at the same T and /i when all 
Thus, the deconfinement phase transition seems to be 



dragged by the chiral phase transition. 
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Fig. 2: The left panel graphs the T dependence of the chiral condensate {qq) / {qq)T=ix=o and the 
Polyakov loop <P near the CEP. The right panel graphs the chiral and Polyakov loop susceptibilities. 




Fig. 3: The left panel graphs the phase diagram on the T—ji plane. The right panel shows positions 

of the CEP in several versions of the NJL model and the PNJL model. The allowed region of the 
/Xg/Te suggested by a lattice simulation is also indicated. The values of Tg and /Xg in several models 
are explicitly given in Table III. 
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It is known [21] that the Polyakov loop susceptibihty has a broad bump in the larger 



T region in addition to a sharp peak in the smaller T region. In Ref. [2l|, the sharp peak 
was interpreted as a reflection of the chiral phase transition and the broad bump was then 
identified with the critical temperature for the deconfinement phase transition. The broad 
bump does not disappear even if the a'^ interaction is introduced. However, we simply define 
the critical temperature with the sharp peak. The critical temperature thus defined for the 
deconfinement phase transition approaches that for the chiral phase transition as fi increases 
to fie, and both the critical temperatures agree with each other when fi > /Xe- This behavior 
is consistent with the generalized Clausius-Clapeyron relation among multiple first-order 
phase transitions [36|. 

The left panel of Fig. 3 shows phase transition lines of six models. Comparing the 
original PNJL and the PNJL+o"^ models, we see that the interaction makes the CEP 
move toward higher T and lower /i. Meanwhile, comparison between the original PNJL 
model and the PNJL model with vector interaction shows that the vector type interaction 
makes the first-order phase transition weak and does the CEP move toward lower T and 
higher fi, that is, in the direction opposite to the case of the interaction. The CEP in the 
PNJL model is always located at T higher than that in the NJL model, even if either the 
vector or the interaction are added to the models. 

The right panel of Fig. 3 shows positions of the CEP in several versions of the NJL and 
the PNJL models. Results of all models are far from the empirical value, but the result 
of the PNJL+o"^ model is closest. It is suggested in the recent lattice analysis js] that the 
possible region where the CEP exists is fie/Te < 2.5. The allowed region is denoted by the 
bidirectional arrow in the right panel. Only the PNJL+cr^ model satisfies this restriction. 
So far we took To=270 MeV, but the pseudo-critical temperature Tc evaluated by adopting 
this value of Tq is somewhat higher than the prediction of a lattice QCD calculation 3?! . 
So, we have rescaled Tq so that the average of the chiral and deconfinement pseudo-critical 
temperatures can agree with that of the lattice QCD at yU = 0. The rescaled Tq is 170 
MeV. The CEP given by the PNJL+a^ model yields /le/Te ~ 1.8 when ro=270 MeV and 
2.1 when To=170 MeV. Both the results satisfy the restriction fie/Tc < 2.5. Moreover, the 
latest lattice analysis points out T^/Tc ~ 0.8. The ratio calculated with the PNJL+a^ 
model is 0.6-0.7 in both cases of To=170 and 270 MeV and consistent with the prediction of 
the lattice analysis jsj. 
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In summary, we have investigated effects of the and the vector interactions on the po- 
sition of the CEP and the interplay between the chiral and deconfinement phase transitions. 
In the case of = 0, the interaction shifts the pseudo-critical temperature of the chiral 
transition to the vicinity of that of the deconfinement transition. As for the CEP, the 
interaction shifts it largely toward higher T and lower /i, while the vector type interaction 
shifts it in the opposite direction. The CEP calculated with the PNJL-|-cr^ model is closest 
to the empirical one and is in good agreement with the restriction /le/Te ^2.5 given by the 
recent lattice analysis. Thus, it is quite interesting to investigate roles of the interaction 
in other thermodynamic quantities such as pressure and quark number density. 
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